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Noncollinear magnetic order is typically characterized by a "tetrad" ground state manifold (GSM) 
of three perpendicular vectors or nematic-directors. We study three types of tetrad orders in two 
spatial dimensions, whose GSMs are SO(3) = 5*^/^2, 5*^/^4, and S'^/Qs, respectively. Qs denotes 
the non-Abelian quaternion group with eight elements. We demonstrate that after quantum dis- 
ordering these three types of tetrad orders, the systems enter fully gapped liquid phases described 
by Z2, Zi, and non-Abelian quaternion gauge field theories, respectively. The latter case realizes 
Kitaev's non-Abelian toric code in terms of a rather simple spin-1 SU{2) quantum magnet. This 
non-Abelian topological phase possesses a 22-fold ground state degeneracy on the torus arising from 
the 22 representations of the Drinfeld double of Qs- 



PACS numbers; 

The search for quantum liquid states has been one of 
the main goals of condensed matter theory for decades. 
There are in general two different routes towards this 
goal, starting from two opposite limits. The first route is 
to start with the quantum limit, say the large- A'^ limit of 
the SU(A) antifcrromagnet, and approach the physical 
system through an 1/N expansion. For instance, the 
1 /N expansion within the slave fermion formalism for the 
SU(iV) antifcrromagnet leads to the valence bond solid 
state with no classical counterpart In our current 
work, we will take a second route towards the liquid state, 
which is by quantum disordering the semiclassical state. 
For instance, it is understood that the valence bond solid 
state (VBS) naturally emerges if quantum fluctuations 
destroy the semiclassical Neel order of a spin- 1/2 system. 
This result is based on the observation that the Skyrmion 
of the Neel order of a spin- 1/2 antiferromagnet always 
carries lattice momentum [2I, [sj . 

In general, the ground state manifold (GSM) of spin 
states can be written as 

GSM = SU(2)/ G, (1) 

where G represents the unbroken subgroup of the SU(2) 
spin symmetry in the ordered phase. G is at least Z2 
for spin-1/2 systems, because physical order parameters 
should be invariant under spin rotation by 2tt. In the 
present paper, we will discuss the quantum disordered 
phases adjacent to semiclassical spin states whose unbro- 
ken symmetry G is a discrete subgroup of SU(2), either 
Abelian or non-Abelian. All these states are "tetrad- 
like" states i.e. the GSM can be represented by three 
perpendicular vectors or nematic-directors (Fig. [IJ . We 
will demonstrate that an exotic non-Abelian topologi- 
cal liquid state can emerge after disordering a tetrad ne- 
matic order of a fairly simple spin-1 system. Non-Abelian 
statistics is a much sought-out phenomenon much dis- 
cussed in particular in fractional quantum Hall systems 
and more recently also in certain topological insula- 
tors (superconductors) In the sequel, we will discuss 
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FIG. 1: Three types of tetrad spin order. Type A phase has 
ground state manifold S0(3), which is equivalent to the con- 
figurations of three perpendicular vectors. Type B phase has 
GSM /Z4, two of the three perpendicular vectors are head- 
less directors. In type C phase all three vectors are headless 
directors. 

in turn three tetrad states, Type A, B, and C. 
- Type A, with G = Z^: Let us first take G = Z^^ and 
thus the GSM is now SU(2)/Z2 = S0(3). S0(3) is pre- 
cisely the tetrad manifold, which corresponds to all the 
configurations of three perpendicular vectors (Fig. [T]A). 
One example of this case is the well understood non- 
collinear spin density wave (SDW), for which the three 
perpendicular vectors A^i, Ni and A^3 that characterize 
the GSM are defined as S{f) = N2 cos{2Q-r)+N3 sin(2Q- 
r); iVi ~ N2 X N3. Here Q is the spij-al wave vector of 
the SDW. It was pointed out in Ref. 6] that if quantum 
fluctuations destroy the noncollinear spin density wave 
(SDW) , one interesting possibility is that the system en- 
ters a Z2 liquid state. On the torus this Z2 liquid ground 
state has a four- fold topological degeneracy Let us 
briefly review why this is the case. The most convenient 
way of parametrizing the manifold S0(3) is by introduc- 
ing CP^ spinor fields z = (^1,22)' as follows: 

Ni - zt^z, N2 RelzHa'^az], N3 - Im[z'icr^CTz]. (2) 

It is straightforward to show that the vectors Na are au- 
tomatically perpendicular to each other after introducing 
the spinor z^. Since all the physical vectors Na are bi- 
linears of Zq, the spinor z^ is effectively coupled to a Z2 
gauge field, which makes z^ equivalent to — Zq. 
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Since the homotopy group 7ri[SO(3)] = Zi, the GSM 
S0(3) supports vortex hke topological defects with a Zi 
conservation law. This type of topological defect is often 
called a vison. Pictorially, the vison can be viewed as 
a configuration in which (for instance) A'^i is uniform in 
space, while and iVa have a vortex (Fig.[2]A). After we 
destroy the ordered state with quantum fluctuations, the 
spinor Za is gapped, but the conservation law of the 
vison still persists. This implies that the disordered phase 
of the noncollinear SDW is equivalent to the deconfined 
phase of Z^ gauge theory, where visons also have a Zi 
conservation law 0, 0] • In this phase the gapped spinor 
Za and the vison have mutual semionic statistics, i.e. the 
wave function picks up a minus sign when Za encircles 
the vison adiabatically. 

As was discussed in Ref. , the transition between the 
ordered phase with GSM SO (3) and the Z2 deconfined 
liquid phase is contiuous and belongs to the 3D 0(4) uni- 
versality class. This is because the bosonic spinor field 
Za can also be viewed as a four component real vector, 
whose order-disorder phase transition belongs to the 3D 
0(4) universality class. The gapped Z2 gauge field does 
not introduce singular corrections in the infrared, i.e. the 
3D 0(4) universality class is unaffected by the the pres- 
ence of the Z2 gauge field Q . Because the physical order 
parameters are bilinears of the spinor z^, they acquire 
a relatively large anomalous dimension as compared to 
the standard order parameters at the Wilson-Fisher fixed 
point of the 0(4) Heisenberg model. Specifically, within 
a five-loop epsilon expansion in d = 4 — e dimensions the 
scaling dimension of these composite bilinears would be 
estimated to be 77^ « 1.37 in (2 -I- 1) dimensions . 
- Type B, with G = Z^: Now let us move to the type- 
B tetrad phase. The GSM can be characterized by 
one vector and two directors, where again all three vec- 
tor/directors are perpendicular to each other (Fig. [l]B). 
It is more convenient to describe this manifold using the 
following slightly different representation of Na'. 

Ni - tr[Z''dZa^],N2 ir[Z^dZa%Ni - tr[ZVZCT^], 

z = (zi, Z2Y = (00 + i(/'3, -</'2 + #i)*- (3) 

Z is a SU(2) matrix, sometimes called the SU(2) slave 
rotor field 23]. Z has an action by SU(2)ioft (left 
multiplication) and by SU(2)i.ight (right multiplication). 
While SU(2)ioft transformations correspond to the phys- 
ical SU(2) spin rotation symmetry, SU(2)i.ight transfor- 
mations contain the gauge symmetry as a subgroup. 

Now let us take iVi a vector, while N2 and are both 
headless directors. In order to make iV2 and headless, 
we can couple Z to a gauge field taking values in a group 
with group elements: 

Z4 = {1, -1, ~ia'}. (4) 
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FIG. 2: The configuration of vison defect in tetrad order A 
and half-vison defect in tetrad order B. 

The gauge field always acts on Z by right multiplication. 
Under the gauge transformation 

Z Z{±i(j''), (5) 

both N2 and reverse direction, while iVi remains in- 
variant. Therefore the type B tetrad phase can be un- 
derstood as the condensate of the SU(2) rotor field Z (or 
spinor z^) when it is coupled to the Z4 gauge field with 
gauge group Z4 from Eq. S) Unlike the type A case, N2 
and are no longer themselves physical order parame- 
ters due to the presence of the Z4 gauge field; rather, the 
physical order parameter Qf = N^N^ - i(M)^ i = 2, 3 
is of quadrupolar type. 

In addition to the vison defect discussed in the type A 
phase, the type B phase also has a "half-vison" defect, 
i.e. the configuration in which A^i is uniform in space, 
while N2 and N3 have a half vortex (Fig. [2]B). This de- 
fect has a logarithmically divergent instead of a confining 
energy because N2 and N3 are nematic directors. We can 
also describe this half vortex as a Z4 gauge flux ±ia^ in 
the condensate of Z. After encircling this flux, Z under- 
goes a gauge transformation as in Eq. [51 and N2 and N3 
reverse their directions. 

When the rotor field Z that couples to the Z4 gauge 
field is gapped out, the system is described by a pure 
Z4-gauge theory - a 'Z4-liquid' phase. If the system in 
this phase is defined on the torus, then there can be four 
different fluxes through each cycle of the torus: 0, 7r/2, 
TT, 37r/2. Each of these different flux combinations cor- 
responds to an independent topological sector. There 
is thus a 16 fold topological degeneracy on the torus. 
Recently, one of the authors of the present paper pro- 
posed that the type B phase is an intermediate phase 
of the Hubbard model on the honeycomb lattice [lo| . 
sandwiched between a fully gapped spin liquid phase and 
a pure Neel order with ground state manifold S'^. In 
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FIG. 3: (a), the topological degeneracy can be counted as the 
number of inequivalent commuting gauge fluxes $i and $2 
through both cycles of the torus. (6) , The ring exchange terms 
Eq. [To] are ring product of gauge field on four links of each 
square, and in the ring exchange terms the links are connected 
in the sequence of the arrows around each plaquette. 



Ref. 10], the vector A^i is the Neel order, whereas the 
directors N2 and Nj, are spin nematic orders. The tran- 
sition between type B tetrad order and the liquid 
phase also belongs to the 3D 0(4) universality class, for 
the same reason as in the type A case. 
- Type C, with G = Qg; Now let us move on to the 
type C phase, whose ground state is characterized by 
three perpendicular nematic directors. This phase can be 
obtained from a system of spin-1 SU(2) quantum spins 
pos sessing both two-spin and four-spin interactions 
11, 12 1 . For a spin-1 system it is often convenient [l3] 
to introduce SU(3) Schwinger bosons bi. When (6*) is 
parallel with (&), the S0(3) spin symmetry is broken, 
while there is no spin polarization on any site, thus the 
system only has nematic order. Since (b*) \\ (b), the 



Schwinger boson b can be rewritten as (bi) = e Ni. Ni 
is in fact a nematic director, because the transformation 
N — J> —N can be cancelled by the transformation 9 6+ 
TT, which is part of the U(l) gauge symmetry associated 
with the Schwinger boson b. The vectors Ni are precisely 
the nematic directors in Fig. [Tp. 

Experimentally it was observed that the triangular lat- 
tice spin-1 material NiGa2S4 has no global spin order 
with time-reversal symmetry breaking at low tempera- 
ture, but it still has gapless excitations with linear dis- 



persion [ij]. It has been proposed[ll|, [l^ that the can- 



didate ground state of this system is characterized by a 
'tetrad' of three perpendicular nematic directors on 
three different sublattices denoted by i — 1,2,3. This 
proposed ground state thus has exactly the same GSM 
as that of Fig. [Tp. The physical order parameter of this 
state is the quadrupolar spin order parameter |24j: 



(6) 



Here Sf is the spin-1 operator on sublattices i = 1,2,3. 
This equation defines 7V°, and it is precisely the ne- 
matic director iV^'^ in the previous paragraph introduced 



through Schwingber bosons [l3|. This GSM is equiva- 
lent to that of the biaxial nematic order 25| of a liquid 
crystal [H, [3| ■ Similar "triatic" nematic spin order was 
also found in numerical work on SU(2) spin-1/2 mod- 
els with both two-spin and four-spin interactions on the 
triangular lattice [17[. 

It is still most convenient to describe this phase with 
the SU(2) rotor variable Z, but now Z is coupled to a 
discrete non-Abelian gauge field taking values in the non- 
Abelian Quaternion group Qs, 



= {±1, ±ia^, ±iay, ±ia''}. 



(7) 



Again, the gauge field acts on the rotor field Z by right 
multiplication. As a consequence of the action of this 
gauge group, Ni in Eq. [3] become headless nematic direc- 
tors. 

Now we will describe this gauge theory based on the 
non-Abelian quaterion group Qs in more detail. Follow- 
ing the general construction in Ref. 0], we define an 8- 
dimensional Hilbert space H on each link (z, ^) of the 
lattice, whose basis elements we denote by Iffi.p). Here i 
denotes a lattice site, fi = x,y a unit vector in a (posi- 
tive) lattice direction, and gi^^ denotes any of the eight 
elements of the group Qs- Now we define, for any group 
element h G Qs, and on every link i,±/i of the lattice, 
operators T"/'^^ and ±^ with the following action on 
the basis vector Ig^,^) residing on that link: 



rp h rri h 



(8) 



and Qj ^ turn out to satisfy the following algebra: 



rph qU 



fh 



Qf rph f^if ^) rphf 



(9) 



The dynamics of the discrete gauge field is given by the 
following 'ring exchange' term 26|: 



ring 



E 

h 



K T: 



X trfG;,., - ,G/,.^-, - 



i+x,y 

G 



■-\-x-\-y ^ — x i-\-y,—y i,x 

^_^Gh^J+H.c. (10) 



Gh is the two dimensional representation Eq. [7] of the 
group element h Cz Qs, and denotes summation over 
all group elements hi^^ S Qs on each link. 

The direction of fi in the ring exchange term on each 
plaquette follows the arrows in Fig. [3l Here we always 
assume K > 0, which favors the gauge flux through each 
plaquette to be 1. 
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The SU(2) rotor field Zi is defined on the vertices i of 
the square lattice. Right- and left- multiplication of Z by 
SU(2) transformations SU(2)right and SU(2)ictt is gener- 
ated by the operators J'^ and J?, satisfying the commu- 
tation relations (see also Ref. II8I ) 

\JkL.A,L\ ^ ^^a^cJ'R^L. \Jr.'A\=^- (11) 

In particular, J£ and act as follows: 



Ze 



Ze- 



J0-Jl 



Ze 



'-Z{\2) 



The quaternion gauge group is a subgroup of the 
SU(2)right transformation. 

The full Hamiltonian with both, rotor and gauge fields 
reads 



ing 



Hr 



Ur ja2 , Ul ja2 

2~ + —-Jl, 



(13) 



This Hamiltonian is subject to the following quaternion 
gauge group constraint: 



f^icr'^ f)^^"^ r^^^ 

^ i,+x^ i,-x^ i. 



(14) 



a = x,y, z. The unitary operator ^'^G^. ^ appearing in 

Eq. [13] is the analogue of the conventional term e*^' f' °' 
where Ai^^ is the gauge potential. The quaternion group 
gauge constraint Eq. [14] generates the following gauge 



transformations on both Z and : 



Zi — )■ ZiGf- 



(15) 



where fi S Qs- The Hamiltonian Eq. [T^] is invariant 
under this gauge transformation. We have formulated 
this model on the square lattice, but generalizations to 
other lattices are straightforward. Again, the quantum 
phase transition between the ordered phase and quater- 
nion liquid phase belongs to the 3D 0(4) universality 
class because the Qs gauge field is always gapped. 

When Ul, Ur > t, the SU(2) rotor field Zi is gapped 
out, and the system is described by the pure quaternion 
group gauge theory Eq. 1101 plus the gauge constraints. 
In the spin Hamiltonian, the rotor field Zi can be gapped 
by turning on the following term on the spin Hamiltonian 
considered in Ref. [Ill Il2| : 



H' = J'^- 



Qj, J'>0. 



(16) 



where Qi is the five-com pon ent quadrupole order param- 
eter introduced in Ref. 11-13]. Eq. is an antifcrro- 



quadrupole interaction between the 2nd neighbor sites 
on the triangular lattice. This term energetically disfa- 
vors the system to form a three sublattice tetrad nematic 



order, and we propose that it will drive the system into 
the phase described by the pure quaternion nonabelian 
gauge theory. This gapped phase is a realization of the 
non-Ahelian toric code phase built on a finite group G, 
proposed by Kitaev [Tj]. In the present case G = Qs- Due 
to the non-Abelian nature of the group Qs, this gauge 
theory is known to possess a rich set of gapped excita- 
tions exhibiting non-Abelian statistics, which are char- 
acterized b y th e representations of the so-called Drinfeld 
double [Ulill of the group Qs- These excitations are 
the following: 

(i) magnetic excitations are located at the centers of 
the plaquettes of the lattice (see Fig. 3), and are charac- 
terized by the product of group elements around a pla- 
quette. Since the product can be taken over different 
closed loops enclosing the same "magnetic flux" , a mag- 
netic excitation is not characterized by a group element g 
itself, but by its conjugacy class Cg = {h^^gh : h e G}. 

(ii) electric charges are located at the vertices of the lat- 
tice (see fig. 3) . An electric charge represents a violation 
of the vertex constraint of Eq. [T3] and corresponds to an 
irreducible representation a of the group G. Transport- 
ing an electric charge a around a magnetic flux Cg along 
a closed path yields the representation matrix D^^^g) of 
the group element g. 

(iii) the most general excitation contains both, mag- 
netic and electric charges (often called a "dyon"), and 
is represented by a pair {Cg,a) as follows: when there 



is no magnetic charge. 



C„=i, then a = a is an 



electric charge, i.e. a representation of the group G. 
However when the magnetic charge associated with a 
"dyon" is not vanishing, i.e. when Cg 7^ Cg=i, its elec- 
tric charge a is an irreducible representation a = a of 
the Normalizer N{g) = {h ^ G : hg = gh} of g (con- 
sisting of all those group elements commuting with g), 
which is in general not the entire group G, but only 
a subgroup thereof. - Let us count the total number 
of excitations for the Drinfeld double of the quarternion 
group Qs- We use the following facts: there are 5 con- 
jugacy classes {+1}, { — 1}, {iic"} where a = x,y, z; the 
number of irreducible representations of any finite group 
equals the number of conjugacy classes; the centralizer of 
any of the three conjugacy classes {izcr"} is the Abelian 
cyclic group .Z4 of four elements generated by icr". Thus, 
there are 5 excitations of the form (Ci, a), 5 of the form 
(C-i,a), and 4 of the form {Cia-<^,a), for each a = x,y,z, 
where d labels the representations of This amounts 
to a total of 5 + 5 + 3x4 = 22 excitations. Since for a gen- 
eral 2D topological field theory the ground state degen- 
eracy on the torus equals the number of topological exci- 
tations ("particles"), this degeneracy is 22 in the present 
case. The fusion rules for these 22 particles can be ob- 
tained from the modular S-matrix (through the Verlinde 
Formula) which, in turn, can be obtained in the standard 
manner [20(1 from the Drinfeld double construction. 
In general, for topological phases which are Drinfeld 
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doubles of a finite group, tlie number of topological sec- 
tors on the torus corresponds precisely to the number of 
commuting pairs of gauge inequivalent magnetic fluxes 
through the two cycles of the torus (Fig. [HJi). (The two 
fluxes need to commute in order to keep the system in 
its ground state.) 

It is important to note that the statistics of the exci- 
tations is only well-defined in the disordered liquid phase 
(described by pure Qs gauge theory). The ordered phase 
with a Z-condensate has gapless Goldstone modes, which 
make adiabatic braiding operations impossible. Another 
important difference between the ordered and the disor- 
dered phase is that these non-Abelian defects have log- 
arithmic divergent energy in the ordered phase, while in 
the disordered phase they all have finite energy. 

Summary: In this work we studied a fully gapped topo- 
logical spin liquid state with non-Abelian excitations. 
Despite its complicated effective model description, we 
propose that such state can be realized by disordering 
a rather simple spin order of a spin-1 quantum SU(2) 
magnet. 
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